Fermi-Liquid Interactions in d-Wave Superconductors 
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This article develops a quantitative quasiparticle model of the low-temperature properties of d-wave 
superconductors which incorporates both Fermi-liquid effects and band-structure effects. The Fermi- 
liquid interaction effects are found to be classifiable into strong and negligible renormalizaton effects, 
for symmetric and antisymmetric combinations of the energies of k t and — k J. quasiparticles, respec- 
tively. A particularly important conclusion is that the leading clean-limit temperature-dependent 
correction to the superfluid density is not renormalized by Fermi-liquid interactions, but is subject 
to a Fermi velocity (or mass) renormalization effect. This leads to difficulties in accounting for 
the penetration depth measurements with physically acceptable parameters, and hence reopens the 
question of the quantitative validity of the quasiparticle picture. 
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There is now considerable experimental evidence that 
the cuprate high T c superconductors exhibit the sim- 
ple power law temperature dependences predicted by the 
quasiparticle picture for their thermodynamic and trans- 
port properties at temperatures well below T c . For exam- 
ple, penetration depth measurements find that the super- 
fluid density exhibits a low-tempepature clean-limit lin- 
ear in n T temperature dependence!!!, in agreement with 
theoryn. The NMR relaxation rate exhibits the ex- 
pected T 3 temperature dependence!]. The predicted ef- 
fect of impurities in giving rise ta a universal thermal 
conductivityoj has been confirmeco. The clean limit spe.- 
cific heat varying as T 2 appears to have been observedD'B. 
Even the electrical transport relaxation rate observed in 
microwave conductivity experimental, which had resisted 
explanation for some time, has-|iiow been explained in 
terms of a quasiparticle pictures. 

Whether or not the magnitudes of the coefficients of 
the above power law temperature dependences are accu- 
rately given by a quasiparticle description is at present 
an open question. A recent study correlating these differ- 
ent coefficients^ concludes that the quasiparticle model 
may be successful here also provided a Fermi-liquid inter- 
action factor multiplying the superfluid density is treated 
as an adjustable parameter. The contention of this article 
is that there is however no Fermi-liquid interaction renor- 
malization of the linear in T contribution to the inverse 
square penetration depth. There is instead a renormal- 
ization by a factor involving the ratio of a band Fermi 
velocity to a Landau quasiparticle Fermi velocity. The 
difficulty now is that a physically unreasonable value of 
this renormalization parameter is obtained from experi- 
ment. This reopens the question of to what extent the 
quasiparticle picture can provide an accurate quantita- 
tive picture of the low temperature behavior of the high 
T c superconductors. Recent debate on correctness of the 
quasiparticle picture, is. also occurring in connection with 
ARPES experimentsHLj, and in connection with the role 
of phase fluctuations of the complex order parameter in 



the determination-of the temperature dependence of the 
superfluid densitycJ. 

The potential importance of Fermi-liquid interactions 
in renormalizing the superfluid density has been empha- 
sized in Refs. |l5| and |l6|. These papers note that Fermi- 
liquid renormalization effects in d-wave superconductors 
can be either strong or weak according as the contribut- 
ing quasiparticles are from the entire Fermi surface, or 
confined the the nodal points where the d-wave gap goes 
to zero. Both expect a strong renormalization effect for 
the superfluid density, whereas this article does not find 
such an effect. 

This article shows that the physics of Fermi liquid ef- 
fects in d-wave superconductors has an interesting sym- 
metry property. This manifests itself when the quasipar- 
ticle energies are separated into parts that are Symmetric 
and Antisymmetric combinations of the energies of the 
+k 1 and — k J, states. (The calligraphic letters S and 
A are used here to emphasize the difference with the 
more usual definition of the symmetric and antisymmet- 
ric combinations with respect to +k j and +k j. states 
common in normal state analyses, e.g. see Eq. 1.32 of 
Ref. H) In the presence of Fermi-liquid interactions, 
the (Symmetric and Antisymmetric corrections to the 
quasiparticle energies obey integral equations that are 
independent of each other, and they are renormalized 
differently. The Symmetric energy corrections exhibit 
strong Fermi-liquid renormalization effects, while the the 
Antisymmetric energy corrections exhibit relatively weak 
temperature-dependent renormalizations that can often 
be neglected. 

Temperature gives a (Symmetric correction to the 
quasiparticle energy because +k f and —k [ states are 
affected in the same way by temperature. A super- 
fluid flow generates an Antisymmetric correction since 
the components of +k and — k along the superfluid ve- 
locity have opposite signs. Also the Zeeman interaction 
generates an Antisymmetric correction because the spin 
| and spin J, contributions to the energy have opposite 
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signs. Thus the superfluid density and the magnetic sus- 
ceptibility are negligibly renormalized by Fermi- liquid in- 
teractions, while the effects of temperature (although rel- 
atively small) are strongly renormalized by Fermi-liquid 
interactions. 

The approach of this article to the inclusion of Fermi- 
liquid interactions in the study of the superconducting 
state follows the intuitively appealing approach of Ref. 
[l|, which is consistent with more formal correlation func- 



tion approaches o£3 Rather than starting with a band 
energy in the absence of electron-electron interactions of 
e\ = h 2 k 2 /(2m) as in Ref. however, this article allows 
e k to be an arbitrary function of k so as to be able to 
account for anisotropic Fermi surface effects. To form a 
Hamiltonian from e\ the substitution Kk — > —iTi\7 — eA/c 
is made. Also, this article studies only equilibrium prop- 
erties, and does not develop a kinetic equation. Other 
studies of Fermi liquid interactions in unconventional su- 
perconductivity include Refs. plp^ . 

The Hamiltonian describing the excitations of the su- 
perconducting state has the following form: 



the addition of the field h£ 



hi 



0, where 

k = de b k /dk, to the Hamiltonian. Herep s = TiV6—eA/c 
is the superfluid momentum, which is assumed to be suf- 
ficiently slowly varying spatially that its gradients can be 
neglected. The velocity v| that appears in h k is the bare 
band velocity, unrenormalized by the electron-electron 
interaction, as noted following Eq. 12 of Ref. and it 
is this same velocity that appears in the expression for 
the quasiparticle contribution to the current density (Eq. 
H below) . On the other hand, the electron-electron inter- 
action contributes to the quasiparticle energy £j. defined 
in and following Eq. |J and hence affects the quasiparti- 
cle velocity vf that occurs in E^jf 1 below. The differences 
in these two velocities have important quantitative con- 
sequences for the interpretation of the penetration depth 
data, as will be seen below. 

The Hamiltonian of Eq. [l] can be used to find the ther- 
mal equilibrium expectation value of the electrical cur- 
rent density operator giving, for the electrical current 



density J = r] g p s 
termined by 



J qp , with the gauge contribution de- 
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-Sef + hf where is the Landau quasipar- 



Here Cfc = & - 

ticle energy relative to the chemical potential (neglecting 
quasiparticle interactions), Afc = 8e k + h kl and is 
the momentum-dependent gap function appropriate for 
rf-wave symmetry. Fermi liquid interactions give a contri- 
bution to the energy of a quasiparticle with momentum 
k, a#4 spin a due to other excited quasiparticles which 
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(c; Cot), and S indicates a varia- 



where a = —a, nka 
tion due to the excitation of other electrons and holes, ei- 
ther by temperature or by the presence of external fields. 
(The factor L~ 2 occurs in Eq. || because the intention is 
to develop a model applicable to superconductivity in a 
two-dimensional copper-oxide plane of a high T c super- 
conductor having area L 2 .) An important step in the 
analysis, as described qualitatively above, is the separa- 
tion of the Fermi liquid interactions into (Symmetric and 
.Antisymmetric parts defined by 



5£k = o [ & feT - fe-fcJ - H = o [ fe fcT + <fe-*i] ■ (3) 



The quantities h k and hf in ft represent generalized 
external fields. For example, in the case of an exter- 
nal magnetic field acting on the orbital motion of the 
electrons, the gap function will acquire a complex phase. 
This phase factor can be removed by a gauge transfor- 
mation, Cka — > CkaGxp{i6), the end result of which is 
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(4) 



and the quasiparticle contribution given by 

ka 

f{Ek,a) is the Fermi-Dirac distribution function, and the 
Ek,a are the Bogoliubov quasiparticle energies defined 
below. The case of an external magnetic field H acting 
on the spin degrees of freedom is described by taking 
h k = ^bH, hf = 0. In both of these cases, the magnetic 
field acts only on the Antisymmetric mode, and has no 
effect on the (Symmetric mode of excitation. 

In addition to causing changes in the energy of a quasi- 
particle (as in Eq. ||), excited quasiparticles can give rise 
to changes in the gap functionLS. There are however pp 
changes that are linear in the superfluid momcntumE3, 
and this effect will therefore be neglected. 

The diagonalization of the Hamiltonian of Eq. |l| gives 

ft = E E k,A,<ff*,<n E ^ = E °* + ^k + h£ k ) (6) 
ka 

where a = ±1, E k = + A 2 , , and the 7^ a are opera- 
tors creating Bogoliubov quasiparticles. 

Later, the energy = \/ £ 2 + A| describing the 

quasiparticle spectrum in the absence of other excited 
quasiparticles is also used. For a d-w&ve superconductor, 
the quasiparticle energy can be parameterizedQ in the 
neighborhood of the Fermi-surface nodal points (see Fig. 

|l|) as E^ = \/{pivf) 2 + {P2V2) 2 , where pi and P2 are 
components of the momentum relative to the nodal point 
in directions perpendicular and parallel to the Fermi line. 
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At low temperatures, only quasiparticles close to these 
four points can be thermally excited. 

Using Eq. || in Eq. ^, keeping only terms up to lin- 
ear order in the fc's and ft.'s, and dropping some terms 
of order ksT / (hkpv 2 ) relative to those kept, yields the 
integral equations 



E, 



(0)' 



and 



Se A - — V f { - 

0e k - L 2 1 kk' 



(-) df 



where &=ffi±ft° w . 
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where the indices 1,2,3 and 4 refer to the four nodes in the 
excitation spectrum, as defined in Fig. |l|. The external 
fields hf at the nodes can be similarly classified. 
The solution of Eq. || now yields 



Se A (T) = -h A F A (T)/[l + F A (T)} 



(10) 



with F A (T) = f A ln(2)k B T/(2Trh 2 v F v 2 ). Here T rep- 
resents any of the irreducible representations present in 
Eqs. ^. The f A, s are defined by 



/a 9 — fit + fl3 + 2fi2 
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FIG. 1. The labelling of the nodes on the Fermi surface 
of YBa 2 Cu306 +a; 

Consider first the (Symmetrical corrections to the 
quasiparticle energies (Eq. |7]) , and assume that there are 
no (Symmetrical external fields other than temperature, 
i.e. hf = (as is the case for the external magnetic 
fields of most interest in this article, which are purely 
Antisymmctrical) . Then the only term driving a nonzero 
contribution to Sef is the term on the right hand side 

proportional to and describing the effect of tem- 

perature. This term is proportional to T 3 , thus giving a 
Sef oc T 3 , and will not be important in contributing to 
the properties of interest at the temperatures satisfying 
fcfiT <C Ao (Ao is the maximum gap). Thus Sef will be 
neglected in calculations below. 

Now from Eq. g which determines the Antisymmetric 
corrections to the quasiparticle energies, it is clear that 
only the values of Set and h A at the Fermi surface nodes 
are relevant to the low energy properties. Also, the so- 
lutions of Eq. U can be classified according to the irre- 
ducible representation of the point group C± v (or 4mm) 
describing a tetragonal copper-oxide plane of a high T c 
superconductor, the independent solutions being 
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are the symmetric and antisymmetric combinations of 
the Fermi-liquid parameters familiar from normal state 
analyses. IHI, and /f 2 for example is f£ k , for k and k' at 
nodes 1 and 2, respectively. 

It is also useful to use Eq. ^ to obtain an idea of how the 
(Symmetrical external fields are renormalized by Fermi- 
liquid interactions. It is clear from Eq. |^ that a knowledge 
of the Fermi-liquid interaction on the entire Fermi surface 
is required and that Sef must be determined on the entire 
Fermi surface. To obtain a rough idea of the nature of 
the solutions, consider a circular Fermi surface of radius 
k F and look for a solution of A g symmetry by considering 

a Fermi liquid interaction ff£; = /Jj , independent of k 
and k' , and a (Symmetrical external field independent 
of k. The solution, which is also independent of k on the 
Fermi surface is 
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1 



F A c 
+ F S A g 



(13) 



where Fj[ 



-- f% kF/{it%VF)- In contrast to the 
Antisymmetrical Fermi liquid parameters F A (T) ob- 
tained above, which go to zero linearly with temper- 
ature in the superconducting state in the clean limit 
(and hence have a dependence on temperature T ex- 
plicitly indicated) , the (Symmetrical Fermi liquid param- 
eter F^ is temperature independent and of approxi- 
mately the same magnitude as the corresponding nor- 
mal state Fermi liquid parameter. The same can be 
seen to be true of the (Symmetrical Fermi liquid param- 
eters corresponding to other irreducible representations 
of Ci v . Note that the ratio of the .Antisymmetrical to 
the (Symmetrical Fermi-liquid F parameters is F A / F s « 
(f A /f s ){k B T/(hk F v 2 )}. 

As noted above, the presence of a superfluid momen- 
tum contributes an _4ntisymmetrical external field to the 
Hamiltonian of Eq. |l]. This external field corresponds to 
the E irreducible representation of C^ v with the p sx and 
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p sy components of n s corresponding to the components 
Ex and Ey of Eq. M. The current density is thus easily 
evaluated using Eq~|| with Eqs. |TJ, [lTj and O. The result 
is J qp = r) qp p s where 



Vq P (T) = - 



2ln2e(v b F ) 2 k B T 
[l + F£{T)]irh 2 v F v 2 



(14) 



Note that the Fermi liquid correction does not alter the 
clean limit linear in T contribution to the r) qp (T), but 
rather makes a T 2 contribution (using (1 + F)^ 1 w (1 — 
F + ...)). Thus there are no Fermi liquid corrections to 
the experimentally measured linear in T contribution to 
inverse square penetration depth. The penetration depth 
A is thus given by 



\- 2 (T) = A~ 2 (0) 



8ln2e 2 2 v F 

— -^-a — k B T ■ 

c 2 h v 2 



(15) 



where a = (vp/vp), and A~ 2 (0) is determined by r\ g 
given in Eq. [|. This has exactly the same form as Eq. 6 
of Ref. [ll], which finds (from a detailed analysis of a num- 
ber of experiments) a 2 — 0.43 for Bi 2 Sr 2 CaCu208 and 
a 2 = 0.46 for YF^CuaOy-^. The conclusions here are 
however completely different from those drawn in Ref. pT| 
which, based on previous theoretical work, considered or 
to be a Fermi liquid correction with a not unreasonable 
value. The conclusion of this article is that the exper- 
imentally determined value of a 2 implies a value of v F 
significantly smaller than vf , which is not physically rea- 
sonable. One of the essential features of strongly corre- 
lated electron systems such as the copper-oxide supercon- 
ductors is that the strong electron-electron correlations 
are expected to produce narrow energy bands and large 
quasiparticle masses, leading to vf <v b F . 

The renormalization of the spin susceptibility due to 
Fermi-liquid interactions can be calculated in a similar 
way. The Zeeman interaction of the spin of an electron 
with the magnetic field contributes an Antisymmetric ex- 
ternal field of A y symmetry to the Hamiltonian. It fol- 
lows that the magnetic moment per unit area of a copper 
oxide plane is 



M = -jl E [/(^m) - m,-i)\ = xh 



(16) 



where 



X(T)= - X °Y^ > Xo(r)= ^ 2 ■ (17) 



l + FHry 



nh vfV2 



Note that here also the low-temperature clean-limit lin- 
ear in T magnetic susceptibility is not changed by Fermi- 
liquid interactions. These affect only terms of order T 2 
and higher in the susceptibility. 

This article has given a detailed description of both 
Fermi-liquid effects and band structure effects within 
the framework of a quasiparticle picture of the low- 
temperature properties of d-wave superconductors. This 



opens the way for a detailed quantitative experimental 
test of the quasiparticle picture of the low-temperature 
properties. A classification of Fermi- liquid effects is 
given that separates the strong renormalization and weak 
renormalization effects according to a symmetry prop- 
erty. The application of the results to the interpretation 
of penetration depth measurements is of particular inter- 
est. This corresponds to the weak (and in fact negligi- 
ble) Fermi-liquid renormalization case, and the ultimate 
conclusion is that the experimental results imply a band 
Fermi velocity which is smaller than the corresponding 
quasiparticle Fermi velocity, which is an unphysical re- 
sult. Clearly there are at present problems with the 
quantitative aspects of the quasiparticle picture of the 
low-temperature properties of the high T c superconduc- 
tors and further study is desirable. 
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